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ABSTRACT 

Mean motion resonances are a common feature of both our own Solar System and of extraso- 
lar planetary systems. Bodies can be trapped in resonance when their orbital semi-major axes 
change, for instance when they migrate through a protoplanetary disc. We use a Hamiltonian 
model to thoroughly investigate the capture behaviour for first and second order resonances. 
Using this method, all resonances of the same order can be described by one equation, with 
applications to specific resonances by appropriate scaling. We focus on the limit where one 
body is a massless test particle and the other a massive planet. We quantify how the the prob- 
ability of capture into a resonance depends on the relative migration rate of the planet and 
particle, and the particle's eccentricity. Resonant capture fails for high migration rates, and 
has decreasing probability for higher eccentricities, although for certain migration rates, cap- 
ture probability peaks at a finite eccentricity. More massive planets can capture particles at 
higher eccentricities and migration rates. We also calculate libration amplitudes and the off- 
set of the libration centres for captured particles, and the change in eccentricity if capture 
does not occur Libration amplitudes are higher for larger initial eccentricity. The model al- 
lows for a complete description of a particle's behaviour as it successively encounters several 
resonances. Data files containing the integration grid output will be available on-line. We dis- 
cuss implications for several scenarios: (i) Planet migration through gas discs trapping other 
planets or planetesimals in resonances: We find that, with classical prescriptions for Type I 
migration, capture into second order resonances is not possible, and lower mass planets or 
those further from the star should trap objects in first-order resonances closer to the planet 
than higher mass planets or those closer to the star. For fast enough migration, a planet can 
trap no objects into its resonances. We suggest that the present libration amplitude of planets 
may be a signature of their eccentricities at the epoch of capture, with high libration ampli- 
tudes suggesting high eccentricity (e.g., HD 128311). (ii) Planet migration through a debris 
disc: We find the resulting dynamical structure depends strongly both on migration rate and 
on planetesimal eccentricity. Translating this to spatial structure, we expect dumpiness to de- 
crease from a significant level at e < 0.01 to non-existent at e > 0.1. (iii) Dust migration 
through PR drag: We predict that Mars should have its own resonant ring of particles captured 
from the zodiacal cloud, and that the capture probability is < 25% that of the Earth, consistent 
with published upper limits for its resonant ring. To summarise, the Hamiltonian model will 
allow quick interpretation of the resonant properties of extrasolar planets and Kuiper Belt Ob- 
jects, and will allow synthetic images of debris disc structures to be quickly generated, which 
will be useful for predicting and interpreting disc images made with ALMA, Darwin/TPF or 
similar missions. 

Key words: celestial mechanics - planets and satellites: dynamical evolution and stability - 
protoplanetary discs - planet-disc interactions - zodiacal dust 



1 INTRODUCTION 

Mean motion resonances (MMRs) occur when two objects' orbital 
periods are close to a ratio of two integers, and a particular combi- 
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nation of orbital angles, the resonant argument, is librating. Exam- 
ples in the Solar System include Neptune and Pluto (3:2 resonance) 
and the inner Galilean moons of Jupiter (4:2: 1 Laplace resonance). 
There are also now numerous examples of suspected or confirmed 
MMRs in extrasolar planetary systems (e.g., GJ 876 b and c in a 
2:1 resonance. iLaughlin & Chambers„2001i) . 
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Mean motion resonances also occur between planets and 
small dust particles, as seen in the Earth's resonant dust 
ring ( lDermottetaDll994i) . Some extrasolar debris discs , such as 
Vega, show evidence of non-axisymmetric clumps (H olland et alj 
Il998l ; lwilneretalj|2002[) . and several authors have attempted to 
model these as arising from a planet's resonant perturbations (e.g., 
iKuchner & Holmanll2003l : IWvatll2003b . 

Although resonant orbits occupy only a small volume of phase 
space, they are common because of a locking mechanism which can 
preserve the resonance once attained. If a particle's or planet's or- 
bital semi-major axis changes due to non-conservative forces, the 
bodies can approach a resonance and then remain trapped there 
even under further action of the non-conservative forces. The asso- 
ciated orbital angular momentum change then drives an eccentric- 
ity change, while the semi-major axis ratio remains approximately 
fixed at the resonance. 

There are many mechanisms by which such a semi-major 
axis change can be dr iven. Early work looked at the tidal evolu- 
tion of satellite orbits ( lGoldreicri[l965h . In a protoplanetary disc, 
planets can mig rate by tidal interaction with the gas disc (see 
IChamberil2009l f or a recent review), an d small planetesimals by 
aerodynamic drag jWeidenschillingl 1 977h . In a gas-depleted debris 
disc, planets can migrate by gravitational sca ttering of planetesi- 
mals ^Fernandez & ldll984l ; lKirsh et alj|2009l) . Interplanetary dust 
drifts towa rds the Sun under the influence of Poynting-Robertson 
(PR) drag jBums et al .11 19791) , an d large bodie s can be moved more 
slowly by the Yarkovsky effect jBottke et ^12006) . At the end of 
a star's main-sequence lifetime, plan etesimals can expe rience aero- 
dynamic drag as the star loses mass toong et alJlioiO ). Moreover, 
for a planet orbiting the secondary component in a binary system, 
formation of a disc following mass transfer from the primary to the 
secondary could trigger renewed planet migration ( .PeretslboiOl) . 

Resona nce capture has be en studied by several authors, go- 
ing back to iGoldreichl ( Il965l) . The regime of adiabatic migra- 
tion, where the migration timescale is much longer than the 
resonant argument's libration timescale, has been studied exten- 
sively analytically using a Ha miltonian model (e.g.. lHenrardlll983 : 
iBorderies & Goldreichlll984l) . With adiabatic migration, capture is 
certain if the particle has an eccentricity below a critical value, and 
probabilistic with decreasing probability as eccentricity increases 
beyon d this. Rapid migration was studied using full N-body mod- 
els by IWvattI ( |2003|, henceforth W03) for the case of a planet mi- 
gr ating into a plan etesimal disc, and using the Hamiltonian model 
bv lOuilIenl (l200a henceforth Q06) for general migration scenar- 
ios. Q06 obtained capture probability as a function of migration 
rate and eccentricity for the Hamiltonian containing a single res- 
onant term, and went on to consider the role played by additional 
resonant terms in affecting capture probability. Such terms can be 
important when the planet is eccentric. 

In this paper we extend this work in a different direction, and 
using the Hamiltonian model with a single resonant term we cal- 
culate capture probabilities, libration amplitudes and offsets for 
particles that are captured, and eccentricity jumps for those that 
pass through the resonance without capture. We validate the model 
against the numerical integrations of W03 and Q06. We also discuss 
the application of the model to various migration scena rios, and dis- 
cuss previous studies in its light, in particular those of iReche et al.l 
( I2OO8 ) wh o investigated c apture of planetesimals by a migrating 
planet, and lPermott et alj ( 1 1994 ) who studied the formation of the 
Earth's resonant dust ring, both using N-body integrations. We 
should like to emphasize the role eccentricity can play in affecting 
capture probabilities and libration amplitudes, which, while long 



understood in the Solar System (e.g.. lMurrav & Dermotl|[l99^ , is 
sometimes neglected in studies of extrasolar planets and discs. 

The Hamiltonian model we use has several advantages over 
N-body simulations: (1) it allows some results to be derived analyt- 
ically; (2) it is faster to integrate numerically than the 3-body prob- 
lem; (3) all resonances of the same order reduce to a Hamiltonian 
of the same form, with fewer free parameters than the three-body 
problem. Once a suite of numerical integrations of the Hamilto- 
nian model is performed, it can be applied to any system, without 
the need for running a different N-body integration every time the 
system parameters are changed. 

The plan of this paper is as follows: In § 2 we describe the di- 
mensionless Hamiltonian model. Readers interested in the details 
may read the Appendix which contains the mathematical deriva- 
tion. In § 3 we summarise how physical parameters relate to the 
dimensionless parameters for test particles. In § 4 we describe the 
results of our numerical integrations. In § 5 we compare the Hamil- 
tonian model to N-body simulations. In § 6 we discuss applications. 
§ 7 summarises our work. 



2 DESCRIPTION OF HAMILTONIAN MODEL 

For most of this paper we consider the circular restricted three body 
problem with a massive planet and a massless test particle orbiting 
a central star with a low mutual inclination. (In the Appendix, we 
derive suitable formulae for the case of two massive planets orbit- 
ing a star, although we do not pursue this further in this paper.) A 
mean motion resonance occurs when the ratio of two bodies' mean 
motions is j : j — k where j and k are integers, and when the asso- 
ciated resonant argument 9 = jX2 — {j ~ k)Xi — kvj is librating. 
Here, A ; are the mean longitudes of the inner and outer bodies, and 
zu is the longitude of pericentre of the test particl^B The integer k 
is the order of the resonance; higher order resonances are weaker 
and often the dynamics are dominated by the low-order resonances. 
In this paper, we consider resonances of first and second order. 

We work with the widely-used Hamiltonian model of mean 
motion resonances (e.g., iMurrav & Dermotllll999l) . This model is 
derived by taking the lowest order term in the disturbing function 
and while suitable for low ecce ntricities other term s may be impor- 
tant at higher eccentricity (e.g.. lLee & Pealell2002 ). Here we sum- 
marise the model and qualitatively describe its behaviour. Readers 
interested in the mathematical derivation are referred to the Ap- 
pendix. In § 3 we summarise the mathematical results. 

2.1 Absence of migration 

First we consider the case with no migration forces. A particle is 
in resonance if its resonant argument is librating, typically about tt, 
rather than circulating. The evolution of the canonical angle k{9 + 
tt) is governed by the Hamiltonian 

H = J'^ +13 J +{-!)'' J''^^ cos ke, (1) 

where J, the generalised momentum conjugate to 6, is proportional 
to the square of the particle's eccentricity. There is one parameter, 
/3, which measures how far the particle is from the nominal reso- 
nance location. All resonances of order k, of any value of j, and 

^ When we refer to inner and outer bodies, their elements are subscripted 
1 and 2 respectively; when we refer to a planet and a test particle, their 
elements are subscripted pi and unsubscripted respectively. 
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whether the particle is interior to or exterior to the planet, can be 
reduced to this form; the relation between the dimensionless vari- 
ables and parameters and physical parameters changes, however, as 
described in the next Section. For first order resonances this is valid 
for coplanar systems and for low inclinations: there are no first- 
order inclination terms in the disturbing function. For second order 
resonances in a non-coplanar system the inclination terms may be 
important, and captu re into a second-order inclination resonance is 
possible (Murray & Dermotl|[l999) . Capture into such resonances 
can be handled using the same formalism but with inclination tak- 
ing the place of eccentricity, and changes to the scaling coefficients 
between the Hamiltonian model and the physical variables. 

Level curves of the Hamiltonian for = 1 (i.e., first-order res- 
onances) are shown in Figure[T] In this plot, the resonant argument 
is measured anticlockwise from the positive a;-axis, and the radius 
is proportional to i.e., to the particle's eccentricity. If /3 is con- 
stant, trajectories are confined to these level curves. For jS > —3, 
there is one fixed point close to the origin, and trajectories circulate 
about this, so no resonant trajectories exist. For jS < —3, there are 
stable fixed points both close to the origin, and removed from the 
origin at 9 = tt. The trajectories librating about this latter point are 
the resonant trajectories. They may have a larger or smaller libra- 
tion amplitude. There are also circulating trajectories at smaller and 
larger eccentricity than the libration region. The curve demarcating 
these regions is the separatrix, shown as a dotted line in Figure[T] 

2.2 Presence of migration 

Now we describe what can happen if there is a non-conservative 
force acting to change the semi-major axis ratio. Migration of either 
the inner or outer body, or both, is mathematically identical: only 
the relative migration rate hi/ a\—a2 / a2 matters. Note that cap ture 
is only possible for converging orbits dMurrav & Dermotll 19991) be - 
cause of the geometry of orbital conjunctions ( Chiang et al.l 2002h . 
so we restrict ourselves to this case. 

Prior to reaching the commensurability, j3 is positive, and all 
trajectories are circulating as shown in the top panel of Figure [T] 
Migration reduces the value of /?, until at /3 = —3 the separatrix 
forms and libration in the resonance becomes possible. As /? con- 
tinues to decrease, the trajectory may either become trapped in the 
resonant region, or continue circulating about the origin. 

In the case of adiabatic migration, progress can be made ana- 
lytically by exploiting the invariance of the area enclosed by the tra- 
jectory, and the criterion for capture is well understood: if the initial 
eccentricity is small, the separatrix forms outside the trajectory, and 
the trajectory is then trapped inside the libration region. Note that 
the trajectories with small J in the top panel of figure[T]deform con- 
tinuously into those in the libration region. If the initial eccentricity 
is large, the separatrix forms inside the trajectory, and initially the 
angle still circulates. However, continued migration causes the area 
enclosed by the separatrix to expand, while that enclosed by the 
trajectory remains constant, so that soon the separatrix hits the tra- 
jectory. Then the area enclosed by the trajectory changes and the 
trajec tory jumps to either the l ibration or the inner circulation re- 
gion dMurrav & Dermotllll99Si) . Although the equations are deter- 
ministic, when one has no information on the phase at which the 
separatrix meets the trajectory, capture can be treated as a proba- 
bilistic event, with the probability of capture decreasing as the parti - 
cle's initial eccentricity is increased dSorderies & Goldreichll984l) . 
If the particle is captured into resonance, it remains so under further 
migration. If the particle is not captured into resonance, circulation 
continues but with a different eccentricity (which is always smaller 
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Figure 1. Level curves of the Hamiltonian for first-order resonances. In the 
absence of migration, trajectories are constrained to lie on these curves. The 
resonant argument is measured clockwise from the origin, and the radial dis- 
tance is proportional to the eccentricity. Top: prior to resonance crossing, 
with /3 positive. All trajectories circulate about the origin. Bottom: after 
resonance crossing, with /3 negative. As well as trajectories where the reso- 
nant argument circulates about the origin, librational solutions are possible; 
these trajectories have been captured into the resonance. The curve demar- 
cating these two regimes is the separatrix (dotted line). The trajectories with 
small J in the upper panel are now those in the libration region. 



for adiabatic migration). If migration continues, however, the par- 
ticle may be captured into another resonance. 

Following capture, the libration continues, with the particle's 
eccentricity increasing as the migration moves the fixed point fur- 
ther from the origin. The centre of libration is also offset slightly 
from the case with no migration. Finally, as we show in the Ap- 
pendix, the libration amplitude decreases slightly, as the angular 
width of the libration region changes while the area remains fixed. 
If the migration force is removed, there will no longer be an offset 
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of the libration centre, and the hbration amphtude will no longer 
vary. 

The case of rapid migration has been less well explored. Q06 
used this Hamiltonian model to derive capture probabilities as a 
function of both migration rate and initial eccentricity. Here we ex- 
tend these investigations to more thoroughly map out the parameter 
space, to find the libration amplitude and offset, and to find eccen- 
tricity changes if capture does not occur. This allows us to build up 
a full picture of what happens as a result of an encounter with a 
resonance. 

Second-order resonances show qualitatively similar behaviour 
to the first-order resonances. 



3 SCALINGS FOR TEST PARTICLES 

Here we summarise how the dimensionless variables of the Hamil- 
tonian model relate to the physical variables and parameters of a 
real star-planet-test particle system. The derivation of these equa- 
tions is given in the Appendix. 



3.1 First-order resonances 

The dimensionless Hamiltonian is 

U = J'^ + pj - J^''^ cose. 



(2) 



For first-order resonances, the dimensionless quantities J and 
P, and the dimensionless time t' , relate to the physical parameters 
as follows: 
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Here, a is the ratio of the planet's and particle's semi-major axes: 
Q = a/spi for an internal particle and a = Opi/a for an external 
particle. The coefficients kj and Ij depend on the particular reso- 
nance and are tabulated in Table lUJ. For test particles outside the 
planet's orbit the "e" columns are to be used; for particles inside 
the planet's orbit the "i" columns are to be used. 

When capture occurs and migration continues, the eccentric- 
ity, libration amplitude and offset evolve with time in the following 
manner: 
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Table 1. Numerical coefficients for conversion between the physical and 
dimensionless variables for first order internal (i) and external (e) test par- 
ticles. Note that the authors believe that any errata in Quillen (2006) have 
been corrected in this paper 
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Table 2. Numerical coefficients for conversion between the physical and 
dimensionless variables for second order internal (i) and external (e) test 
particles. 



Eccentricity is pumped up, while the libration amplitude and offset 
are reduced. Note that the value of the libration amplitude is not 
determined analytically, but its dependence on time can be. 

3.2 Second-order resonances 

For second-order resonances, we write the dimensionless Hamilto- 
nian 



.r + pj + jcos2e. 

The scalings for J, /3 and t' are then as follows: 
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The coefficients are tabulated in Table l|2j. 

When capture occurs and migration continues, the eccentric- 
ity evolves according to Equation |6l libration amplitude and offset 
evolve with time in the following manner: 
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This is similar to first-order resonances, but the time dependence 
of the libration amplitude and offset is stronger. Again, the value 
of the libration amplitude is not determined analytically, but its de- 
pendence on time can be. 



4 NUMERICAL INTEGRATION 
4.1 Capture probabilities 

4.1.1 First order resonances 

For our numerical investigations, we show plots of capture prob- 
ability, libration amplitude and offset, and eccentricity jump, in a 
region log |/3| G [—3, 1), log Jo G [—2, 2). Here Jo is the initial 
value of the momentum J. Note that higher-order terms in the dis- 
turbing function will become important for e > 0.1, i.e., Jq > 1 or 
100 for a Jupiter-mass or an Earth-mass planet respectively. Note 
also that an eccentricity of 1 is attained at Jo ~ 100 or 10'*, for 
Jupiter-mass and Earth-mass planets respectively. Hence, results on 
the extreme right-hand side of the following figures may not be ac- 
curate for massive planets. 

We consider the Hamiltonian given by Equation (|2}. For the 
integration we transformed to Poincare's canonical Cartesian vari- 
ables (equation lA45b in order to avoid a singularity in the equations 
of motion for small J. The Hamiltonian is then 



H 



(15) 



The equations of motion arising from this Hamiltonian (i = 
—dH/dy, y — dT-L/dx) were integrated numerically with a 
Runge-Kutta routine with Cash-Karp spacings and adaptive step 
size (iPress et al...2007') . We also vary /? by including the equation 
$ = constant to impose migration; this changes semi-major axis 
without affecting eccentricity. Note that it is sufficient to consider 
linear variation in /3 since it is only the migration rate at the in- 
stant of resonant passage that affects capture probability, so long 
as the migration is sufficiently smooth. However, if the migration 
rate changes rapidly, such as if the planet is experiencing turbu- 
lent torques, then the model does not apply. The initial value of 
/3 was chosen to be 15 ± 1 with a uniform random distribution. 
It was found that when all particles were started from exactly the 
same value of /3 the capture probability varied in a very rapid man- 
ner with Jo 01 $. This appears to be related to the phase at which 
particles enter the resonance; the distribution of phases does not re- 
main uniform as the system evolves away from the initial uniform 
distribution. While perhaps mathematically interesting, this is not 
physically useful, so we randomised over this parameter. 

We varied the initial momentum Jo and the migration rate j3. 
For each point in parameter space we integrated 100 trajectories, 
with the resonant argument chosen from a uniform distribution over 
[0, 2n) and initial /? as described above. A trajectory was classed 
as a capture if all of 1000 output values of 6 were within 47r/5 ra- 
dians of the centre of libration. This avoids any misclassification: 
choosing a larger libration amplitude would risk misclassifying cir- 



culating trajectories as captures. The top panel of Figure ([2} shows 
capture probabilities as a function of Jq and 43- 

In Figure (2] top panel, we see that capture into resonance is 
guaranteed for small initial eccentricities and migration rates. For 
low migration rates we are in the well-studied adiabatic regime. 
For low eccentricities capture is certain since the separatrix forms 
around the initial trajectory. For high eccentricities the separatrix 
forms inside the initial orbit and expands to meet it as the migra- 
tion continues. Capture then is probabilistic, with a probabili ty that 
decreases as the initial eccentricity increases ( lHenrard[T982h . 

For low eccentricities, with Jo < 1.3, capture is certain if the 
migration rate is low and impossible if it is high. The width of the 
transition region from capture to no capture increases with eccen- 
tricity. In the limit of low eccentricity, the transition occurs at a crit- 
ical migration rate of |/3| ~ 2.1. Certain capture occurs for low mi- 
gration rates up to Jo ~ 1.3. For higher eccentricities (Jo > 1.3), 
capture is always probabilistic, with a capture probability that is not 
strongly dependent of migration rate; however, if migration rate is 
too high, then capture is still impossible. Interestingly, the maxi- 
mum migration rate allowing capture is slightly higher for higher 
eccentricities. This is because this maximum migration rate is gov- 
erned by the resonant libration period, which decreases with higher 
J (Eq. |A49l >. For a given migration rate above the critical rate, the 
capture probability peaks at a finite eccentricity. 

Vertical line cuts through Figure (2] bottom panel, give the 
same results as shown in figure 2 of Q06, who integrated the same 
Hamiltonian with a different integrator, providing a check on our 
numerical implementation. 

In subsequent sections we shall sometimes be concerned with 
particles that have dimensionless momenta Jo > 100, for example 
when considering planets of around Earth mass or lower. To calcu- 
late the capture probabilities at these high momenta, we extrapo- 
lated the results plotted in Figure [2] At high momentum (Jo > 10) 
the capture probability is independent of migration rate for most 
of the range of migration rates considered, and a regression fit af- 
ter poohng the data for Jo G [3.2, 100) and $ G [0.001, 4.2] 
gives p — 0.499JJ''''^^^ for the capture probability. This agrees 
with the analytical re sults for the adiabatic case where p oc e~^^'^ 
( lDermottetalJl988l) . 



4.1.2 Second-order resonances 

We used Poincare variables and same integrator as for the first- 
order resonance to integrate the second-order Hamiltonian 



l)x^ ^ il3~l)y' 



(16) 



Again we randomised the initial angle and distance to the reso- 
nance. The integrations were performed for log \$\ G [—1.5, 1.5), 
log Jo =G [—9, 1.6875). Higher order terms will become impor- 
tant for J of several hundred. The capture probabilities are shown 
in Figure bottom panel. 

The overall picture is similar to that for first order resonances: 
Certain capture at low migration rate and eccentricity, no cap- 
ture at high migration rate and low eccentricity, and probabilistic 



^ The data used to create figure |2] and the other other 
contour plots in this paper are available on-line at 
|http : / / www . ast . cam . ac . uk / - a jm233/ and the journal website 
and may be used provided that this work is cited. 
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Figure 2. Top: Capture probabilities for the first-order resonance, as a 
function of rescaled eccentricity (x-axis) and migration rate (y-axis). Con- 
tours are at 1%, 25%, 50%, 75% and 99% probability. Also shown (green 
crosses) are the dimensionless migration rate and eccentricity for a par- 
ticle with Keplerian eccentricity 0.01, migrating at lAUMyr^^ into ex- 
terior resonances with an Earth-mass planet orbiting a Solar-mass star at 
1 AU. The 2:1 and 3:2 resonances are labelled. Note that the axis of or- 
dinates is condensed at the bottom of the plot. The bottom right comer 
(log Jo G [1.75, 2), log $ £ [—3, —1)) was not integrated due to long in- 
tegration times. Bottom: Capture probabilities for second-order resonances. 
Contours are at 1%, 25%, 50%, 75% and 99% probability. 



capture at low migration rate and high eccentricity. In the adia- 
batic hmit, the transition from certain to probabilistic capture oc- 
curs at Jo ~ 0.8. At high eccentricities, capture probabihty does 
not depend strongly on migration rate for the rates considered, 
and we find p = .867Jrr°'^^, agreeing with analytical studies 
jPermott et alj|l98^ which give p oc e^^. In the limit of low ec- 
centricity, the transition from certain to impossible capture does 



not appear to converge to any value of /3, unlike for first order res- 
onances, being quite strongly dependent on the migration rate, par- 
ticularly for Jo > 0.001. At higher eccentricities, the maximum 
migration rate allowing capture increases more rapidly with J than 
for first order resonances. This is due to the stronger dependence of 
the libration period on J (Eq. |A52t . 

Vertical line cuts through Figure (2] bottom panel, give the 
same results as shown in figure 3 of Q06. 



4.2 Libration amplitude and evolution in resonance 

The analytical results in section |A3] show that libration amplitude 
and offset decrease with time, so we need to specify a time at which 
to measure the libration amplitude and offset. First we check the ac- 
curacy of these analytical results. For first order resonances, the top 
panels of Figure[3]and Figure|4]show the evolution of libration am- 
plitude and offset respectively, for one trajectory from each corner 
of parameter space (Jo = 0.1, /3 = 0.1; Jo = 0.1, /3 = 2.5; 
Jo = 10, /3 = 0.1; and Jq = 10, /3 = 7.5). The points are the 
numerically determined amplitudes and offsets from an extended 
integration of these trajectories. For the offsets, the lines are the 
analytical solution given by Equation lA44l For the amplitudes, the 
lines have a slope of — 1/8 as determined analytically, and a nor- 
malisation given by the relevant shorter integration from the grid 
described above. Agreement is generally good, for both fast and 
slow migration, even though the result was derived assuming adi- 
abatic conservation of area. This demonstrates the accuracy of the 
analytical solution, and also that the output from the grid of inte- 
grations is sufficient to correctly describe the long-term behaviour. 

The bottom panels of Figures [3] and |4] show the analogous 
plots for second order resonances (trajectories taken with Jo — 
0.032,/3 = 0.032; Jo = 0.032, $ = 2.4; Jo = 13, /3 = 0.032; 
and Jo = 13, P — 21). Again we see good agreement between the 
analytical results and the numerical integrations. 

The libration amplitudes and offsets for first and second order 
resonances for a range of Jq and /? are shown in figures |5] and [6] 
Note that all amplitudes and offsets are scaled to what they would 
be at /3 = — 100 using the time dependence given in equations |A44] 
|A48l|A50l and |A51| Large amplitude librations can occur as a re- 
sult of either large initial momentum or, to a lesser extent, faster mi- 
gration. For first order resonances, at low momentum the libration 
offset is not dependent on Jo and increases with m igration rate, as 
predicted by analytic theory jPermott et alj[l994h . At higher mo- 
menta the offset decreases with increasing Jo, although here the 
data are noisy due to low numbers of particles being captured. For 
second order resonances the offset increases with increasing /3, and 
as Jo increases the offset decreases, passing through a minimum 
before increasing again at high Jo. 

Note that the offset is only present during migration: if the 
migration force is removed, such as if a perturbing gas disc has 
dissipated, then libration will be about tt exactly with no offset. 



4.3 Failure to capture 

In the event that the particle not be trapped in resonance, its ec- 
centricity changes as it passes through the resonance. In the case of 
adiabatic migration, the eccentricity is reduced ( Murray & Dermot^ 
On the other hand, we find that when the particle fails to cap- 
ture due to fast migration rate, the eccentricity can change substan- 
tially in either direction. The top panel of Figure |7] illustrates how 
the eccentricities are changed when passing through a first-order 
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Figure 3. Top: Evolution of libration amplitudes (in degrees) with time, 
for the first-order resonances. Lines show the analytical result from § I A3 1 
and points show the numerical results. Bottom: The same, for second-order 
resonances. 

resonance. The colour scale and solid contours show the final ec- 
centricity as a fraction of the initial eccentricity. We see agreement 
with adiabatic theory in the case of slow migration: all particles' 
eccentricities are driven down. In Figure |7] we also show contours 
of (mean change in momentum)/( standard deviation of change in 
momentum) to illustrate how broad the distribution of eccentricity 
changes is relative to the mean. For slow migration rates the dis- 
tribution is very narrow, and all particles behave in the same man- 
ner. In contrast, for faster migrations and high momentum we see 
that not only can the average change in eccentricity be positive, but 
that the distribution of changes is relatively broad, so that many 
particles lose eccentricity despite the mean being an increase. In 
this case, the behaviour is highly stochastic. For fast migration and 
low momentum, the particles' eccentricities are all driven up. We 
note furthermore that the actual distribution of eccentricity jumps 
is highly non-Gaussian, being either unimodal and strongly skewed 



Figure 4. Top: Evolution of the offset of libration centres (in degrees) with 
time, for the first-order resonances. Bottom: The same for second-order 
resonances 

or even bimodal. Note that for low initial eccentricities the mean 
eccentricity jump decreases as the migration rate is increased, pre- 
sumably because the resonance has less time to affect the particle's 
orbit before being crossed. 

The bottom panel of Figure |7] shows an analogous plot for 
the second-order resonances. The behaviour is qualitatively simi- 
lar: with slow migration the eccentricity of all particles is decreased 
on passing through the resonance, but for fast migration and low 
eccentricity the eccentricity is pumped up, while at higher eccen- 
tricity the eccentricity may jump either up or down. 



5 MODEL VALIDATION 

We compared line cuts through our figures of capture probabil- 
ity to figures 2 and 3 of Q06, demonstrating that our integrations 
of the Hamiltonian are accurate. To demonstrate that the Hamil- 
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Figure 5. Top: Mean libration amplitudes for the first-order resonances, in 
degrees. Amplitude is that at /3 = -100. Contours are at 15°, 30°, 60°, 
and 120°. Note that the axis of ordinates is condensed at the bottom of 
the plot. Bottom: Mean libration amplitudes for the second-order reso- 
nances, in degrees. Amplitude is that at /3 = —100. Contours are at 
7.5°,15°,30°,60°,andl20°. 



Figure 6. Top: Mean libration offsets for the first-order resonances, in de- 
grees. Offset is that at /3 = -100. Contours are at 4°, 8°, 16°, and 24°. 
The data are noisy at high Jq where there are few particles captured. Note 
that the axis of ordinates is condensed at the bottom of the plot. Bottom: 
Mean libration offsets for the second-order resonances, in degrees. Offset is 
that at fl = -100. Contours are at 0.25°, 0.5°, 1°, 2°, and 4°. 



Ionian model itself coiTectly describes the behaviour of the three 
body problem, we now compare our results with the N-body sim- 
ulations of W03, who considered the case where a massive planet 
was migrating outward into a disc of test particles. We conduct a 
Monte-Carlo simulation, creating a population of particles with the 
same distribution as in the N-body integration and converting their 
parameters into the dimensionless ones to determine the outcome 
of the resonance passage process. This process is very quick: each 
point in Figure [8] represents 100 particles, so there are around 10® 
particles in total, and the process took only a few seconds. Hav- 
ing generated our samples we plot their capture probabilities and 



amplitudes and offsets against migration rate, and also show the 
empirical fitting formulae from Wyatt's N-body simulations. 

The top panel of Figure [8] shows the capture probability as a 
function of planetary migration rate for the 3:2 external resonance. 
Test particles were located at 60 AU, with eccentricities in the range 
[0,0.01]. Hence the population for each point in (d, m-pi) space is 
taken form a horizontal cut through the dimensionless parameter 
space, up to a certain maximum value of J. The stellar mass is 
2.5 Solar masses and the planet mass ranges from 1 to 300 Earth 
masses. The probabilities obtained from the Hamiltonian model 
show excellent agreement with those obtained from N-body inte- 
grations (see figure 1 of W03); in particular we reproduce the in- 
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Figure 7. Top: Change in eccentricity when passing through a first-order 
resonance, in the event of a non-capture. Colour gradient and associated 
solid contours show the mean final momentum J as a fraction of the initial 
Jo, so that particles on the 1.00 contour experience no mean change, par- 
ticles with ratio greater than 1 .00 are pumped up on average, and particles 
with ratio less than 1.00 are cooled on average. The dotted contours show 
the mean change divided by the standard deviation. At low migration rates 
the distribution is narrow and all particles behave in essentially the same 
way — they are cooled — but at high migration rates there is a large spread in 
eccentricity jumps. Bottom: The same, for second order resonances. 



creasing sharpness of the transition from certain capture to certain 
failure to capture as the planet mass is increased. This is because for 
low mass planets, an eccentricity of only 0.01 is sufficiently non- 
zero to affect the trapping probability, putting the particles in the 
region of parameter space where the transition from certain capture 
to impossible capture as migration rate increases begins to broaden. 
Models involving capture into resonance with low mass planets will 
therefore need to pay careful attention to the initial conditions of the 
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Figure 8. Top: Capture probabilities as a function of planetary migration 
rate for the 3:2 external resonance. Solid lines are the fitting formulae 
from the N-body simulations of W03, and points are from the Hamilto- 
nian model. From left to right, planet masses are 1, 3, 10, 30, 100 and 300 
times Earth's mass. Bottom: Capture probabilities as a function of plane- 
tary migration rate for the 5:3 external resonance. Solid lines are the fitting 
formulae from W03, and points are from the Hamiltonian model. From left 
to light, planet masses are 30, 100, 300 and 1000 Earth's mass. 



planetesimal disc, since the distribution of eccentricities, even in a 
dynamically cold disc, can affect resonant trapping behaviour. In 
contrast, capture into resonance with high mass planets will follow 
the same behaviour for all small eccentricities. 

The capture probabilities for second-order resonances, too, 
show excellent agreement with N-body simulations. Compare the 
bottom panel of Figure[8]with figure 4c of W03, for the case of the 
5:3 resonance. The second order resonances are much weaker than 
the first order ones, and we cover the same range of migration rates 
with more massive planets: from 30 to 1000 earth masses. 

In the top panels of figures |9] and [TO] we compare the off- 
sets and amplitudes from the Hamiltonian model with the results 
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Figure 9. Top: Libration offsets for the 3:2 external resonance. Solid lines 
are the fitting formulae from W03, and points are from the Hamiltonian 
model. From left to right, planet masses are 1, 3, 10, 30, 100 and 300 
times Earth's mass. Bottom: Libration offsets for the 5:3 external reso- 
nance. Solid lines are the fitting formulae from W03, and points are from 
the Hamiltonian model. From left to right, planet masses are 30, 100, 300 
and 1000 times Earth's mass. 

from IWvat for the 3:2 external resonance. We see quali- 

tative agreement between the Hamiltonian model and the N-body 
fits, although quantitative agreement is not exact. We see the li- 
bration offset steadily increasing with migration rate, although that 
from the model is somewhat lower. The libration amplitudes are, 
for planets more massive than 1 M®, constant for small migration 
rates before increasing at higher migration rates, although for Earth 
mass planets the amplitude is independent of migration rate. At a 
given migration rate, capture by a less massive planet will result in 
a larger libration amplitude. The fitting formulae from W03 shown 
on this plot were found by aggregating over all planet masses, so the 
Hamiltonian model provides a more accurate prediction for a given 
mass. The difference between the Hamiltonian and N-body results 
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Figure 10. Top: Libration amplitudes for the 3:2 external resonance. Solid 
lines are the fitting formulae from W03, and points are from the Hamilto- 
nian model. From left to right, planet masses are 1, 3, 10, 30, 100 and 300 
times Earth's mass. Bottom: Libration amplitudes for the 5:3 external res- 
onance. Solid lines are the fitting formulae from W03, and points are from 
the Hamiltonian model. From left to right, planet masses are 30, 100, 300 
and 1000 times Earth's mass. 

may be due to W03 not using a consistent duration of migration 
post capture; we have taken the migration to continue for 0.5 AU 
in all cases. We also have a semi-analytical motivation for Wy- 
att's finding that the libration amplitude depends on dm~^ '^: the 
libration amplitude depends on the dimensionless migration rate $, 
which is in turn proportional to am^^^^ (equation [3}. 

The bottom panels of Figures [9] and [TO] show the analogous 
plots for the 5:3 resonance. Again the model reproduces the qual- 
itative behaviour of the N-body results, with the model this time 
yielding slightly higher offsets, and similar differences between the 
model and N-body as in the case of first order resonances. Again, 
higher libration amplitudes result from lower planet masses at a 
given migration rate. 
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6 DISCUSSION 



therefore have 



6.1 General comments 

To summarize the previous sections, we have shown that capture 
into a resonance is certain at low eccentricities and migration rates, 
impossible at fast migration rates and low eccentricities, and possi- 
ble but not certain for high eccentricities. The critical eccentricity 
at which capture becomes probabilistic is lower for lower planet 
masses: around 0. 1 for a Jupiter mass planet and around 0.01 for an 
earth mass planet, for first order resonances. Furthermore, at lower 
eccentricity, the width of the transition from capture to no capture 
with increasing a is also affected. The effects of this can be seen in 
Figure[8] Eccentricities of order 0. 1 may not be uncommon in plan- 
etesimal belts: bodies in the cl assical Kuiper Belt c urrently have a 
mean eccentricity of about 0. 1 jLuu & Jewitll20o3) . and eccentric- 
ities in planetesimal belts can be pumped up above 0. 1 when the 
largest bodi es grow significantly and gravitationally perturb their 
neighbours jKenvon Sc Bromlevll2008[) . Thus, the role of eccentric- 
ity cannot be neglected when investigating resonant encounters. If 
capture into a resonance occurs, the resulting libration amplitude 
is larger if initial eccentricity was larger, and, for small initial ec- 
centricities, also increases with increasing migration rate. Libration 
centre is offset by an amount dependent only on migration rate, in 
agreement with analytic theory, unless the dimensionless momen- 
tum is very high. If a particle fails to be captured, the eccentricity 
is invariably driven down for slow migration rates, but can jump up 
or down for high migration rates. 

Before discussing specific scenarios, we make some general 
comments. 

For any migration scenario, if migration lasts sufficiently long, 
a succession of resonances will be encountered, if capture fails 
to occur at the first resonance. For low eccentricities, the criti- 
cal migration rate for capture increases with j, so a particle fail- 
ing to capture into the 2:1 resonance may be captured into the 
3:2, 4:3, etc. Indeed, the critical migration rate increases with- 
out bound as j — >■ oo, so it would appear that ultimately a par- 
ticle will be captured into a resonance very close to the planet. 
However, such resonances are not stable. Close to the planet the 
resonances overlap, forming an unstab le region of c haotic be- 
haviour of width ^ 1.4api(mpi/m^)2/'' l lWisdonJl98"oh . Thus the 
closest first-order resonance which is stable is given by jmax = 
18(mpi/m0)~^''^(m*/m0)^''^. Even if particles were captured 
into these resonances, they would be quickly removed. Hence for a 
given planet mass there is a maximum j which will trap particles, 
and conversely, for a given resonance there is a maximum planet 
mass, which is determined by the width of the chaotic zone. This 
is illustrated in Figure [TT] The top and bottom panels are for cap- 
ture of particles exterior to and interior to the planet, respectively. 
The diagonal lines show the critical migration rate for capture into 
resonances of a given j as a function of planet mass. This figure 
generalises figure 5 of W03 to resonances of higher j, illustrating 
the power of the Hamiltonian model. 

The lines in FigurefTTIare terminated by the chaotic zone pre- 
scription described above. For example, a 1 planet at 1 AU 
will be able to capture exterior low eccentricity particles in res- 
onances up to the 18:17. To calculate the upper envelope, we 
know that the critical dimensionless migration rate is /3 ~ 2.1 
(^ 14.1.11 ). From table [T] we find a fit Ij « 5.4j"^-*, so the rel- 
evant coefficient for the innermost stable resonance is Ij^i-m,^ ~ 
0.0297(mpi/m®)"-^"(m*/mQ)-°-^^'*. Using Equation H] we 
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above which capture is impossible into any resonance, unless the 
eccentricity is sufficiently high. This equation holds for both exter- 
nal and internal particles. 

While Figure[TT]is plotted for a planet at 1 AU, similar results 
can be obtained for other distances: the critical migration rate for 
a given resonance decreases (oc l/yTipT, from equation[3j and the 
critical migration rate for capture into any stable resonance has the 
same dependence. Thus, slightly lower migration rates are required 
for planets further from the star to capture particles than planets of 
the same mass in smaller orbits. The effect of this will depend on 
the migration mechanism: for Poynting-Robertson drag, the migra- 
tion rate is inversely proportional to a, so lower j resonances will 
capture particles when the planet is further from the star. On the 
other hand, for Type I migration in a disc with a MMSN power law 
index, the migration rate is independent of a, so capture will be into 
higher j resonances if the planet is further from the star. 

We have ignored the effects of eccentricity in the previous 
paragraph. This plays an ambiguous role: at low migration rates 
eccentricity reduces the probability of capture. For first-order res- 
onances, the critical eccentricity for the transition from certain 
to probabilistic capture decreases with increasing j (Figure I12t . 
For second-order resonances the critical eccentricity does not de- 
pend strongly on j. However, at high migration rates, higher ec- 
centricity can make capture possible where capture was impossi- 
ble at low migration rates; this effect is particularly noticeable for 
second-order resonances. Furthermore, resonances can cause par- 
ticles' eccentricities to change as the particles pass through res- 
onance and fail to be captured , and the cause of migration itself 
often damps eccentricities too teums et alJl 19791 : iTanaka & Ward 
|2004[ although disc-planet tor ques may also pump eccentricity, 
e.g.. lMoorhead & Adamsll2008l) . 

If a particle encounters a succession of resonances, then, ig- 
noring changes in eccentricity due to resonant encounters or other 
effects, the dimensionless migration rate will decrease and the di- 
mensionless momentum increase at each subsequent resonance, 
simply from the changing scaling coefficients due to the changing j 
of the resonance (Equations[3]l. We indicate such a locus of (Jo, j3) 
on Figure |2] as green crosses, for a particle with e = 0.01, a = 
1 AU Myr^ ^ for particles encountering external resonances with an 
Earth mass planet orbiting a Solar mass star at 1 AU. Here the par- 
ticle is migrating too fast to be captured into the 2:1 resonance, 
but is very likely to be captured into the 3:2, 4:3 and subsequent 
resonances. Capture probability decreases for higher j. 

Our treatment does not include any eccentricity belonging 
to the massive planet. Correctly accounting for this requires in- 
cluding in the Hamiltonian another resonant term for the reso- 
nanc e involving the p l anet's longitude of pericentre (Q06). How- 
ever, |^uiiien&F^ ( l2006h showed that the strength of this reso- 
nant term can be reduced if the particle's osculating Keplerian ec- 
centricity is replaced by the free secular eccentricity, which is ap- 
proximately the same as the planet's eccentricity for particles close 
to the planet. Hence, for a first approximation, in the subsequent 
discussion we consider the impact of the planet's eccentricity as 
though it were the planetesimals'. 
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Figure 11. Top: Critical migration rate at 1 AU for capture of low- 
eccentricity particles drifting into resonances exterior to a planet as a func- 
tion of planet mass. Stellar mass is Solar. At migration rates faster than this, 
capture is impossible at low eccentricity. The critical migration rates for 
first order resonances are shown as solid lines. The 2:1 and 3:2 resonances 
are labelled; the higher j resonances move away from these monotonically. 
The line for a given resonance is terminated when that resonance becomes 
unstable according to the resonance overlap criterion. We also show second- 
order resonances as dotted lines; only up to jr = 21 are shown for clarity 
(13 = 0.3 was taken as the critical migration rate, although this depends on 
eccentricity; see Figure|2] $ = 0.3 is the maximum of the envelope of the 
100% capture region.). The right-hand axis of ordinates shows the equiva- 
lent dust grain size for migration under PR drag (see ^ 16.3.21 . Bottom: The 
same, but for trapping of particles interior to a migrating planet. We show 
typical migration rates for planets embedded in gas discs at 1 AU (dashed 
lines). Mass-dependent Type I migration occurs for low mass planets and 
mass-independent Type II migration for high mass planets; see ^ 16.21 



Figure 12. Critical eccentricity for certain capture of slowly migrating par- 
ticles into external resonances as a function of planet mass. The star is of 
Solar mass and the planet is located at 1 AU. At eccentricities greater than 
this capture is possible but becomes less probable with increasing e. Criti- 
cal eccentricities for first order resonances are shown with solid hnes. The 
2:1, 3:2, and 4:3 resonances are labelled; the higher j resonances move 
away from the 3:2 and 4:3 monotonically. The line for a given resonance 
is terminated when that resonance becomes unstable according to the reso- 
nance overlap criterion. Critical eccentricities for second order resonances 
are shown as dotted lines. 



6.2 Formation of resonances during planet formation 

Here we discuss the capture of planets into resonances due to mi- 
gration through a gas disc. Due to the complexity of the problem, 
this Hamiltonian model can only give us a crude understanding of 
this scenario. We have neglected the mass of one planet, so when 
the planets are of comparable masses it may be necessary to in- 
tegrate the equations of motion for two planets (Equations IA16I 
and lA28t . In spite of this limitation, the Hamiltonian model may 
give us some basic insight into the capture process. 

There are several scenarios which may arise in a protoplan- 
etary disc. For example, a giant planet undergoing type II mi- 
gration may ca pture planetesimals an d low-mass embryos inte- 
rior to it (e.g.. iFogg & NelsonI 120071) . The outer planets of our 
o wn Solar system may ha ve once been in a chain of resonances, 
as iMorbidelli et al.l ( |2007|) suggested for initial conditions for the 
Nice m odel of the outer Solar System, and which lPierens & Nelsod 
( I2OO8I) suggest may be a general outcome of migration of Jupiter 
and Saturn mass planets. Lower-mass super-Earths migrating in- 
wards can encounter th e ex t erior resonances of a giant planet (e.g., 
Podlewska & Szuszki ewiczi 12009) . On the other hand, there may 
be a tendency for multi-plan et systems to have the mo re massive 
planet on an external orbit jMorbidelli & Cridall2007l) . Here we 
briefly discuss this latter case, where an inwardly migrating planet 
is capturing smaller planets or embryos in interior resonances. 

On Figure [TT] (bottom panel) we have shown the typical mi- 
gration rates due to Type I migration (which is mass-dependent) 
and Type II migration (mass-independent) for typical disc proper- 
ties (MMSN, disc viscosity parametrised by a = 0.01, scale height 
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h = 0.025) at 1 AU ( IChambers|[2009lPi We see that, for the 2:1 
resonance, Type I migration at this location in this disc is always 
too rapid allow capture. Capture into the 2:1 resonance could occur 
however if the planets were located closer to the star at the epoch of 
capture. Although the Type I migration rate decreases with planet 
mass, the critical migration rate decreases more strongly, so cap- 
ture into this resonance is not possible with these disc parameters. 
Only more massive planets moving slowly under Type II migration 
can capture into the 2:1 resonance, which will be an unavoidable 
outcome of migration in this disc if eccentricities are low, unless 
planets formed closer to each other than the 2: 1 resonance. Lower 
mass planets can only capture into higher j resonances. Indeed, we 
should see lower mass planets in higher j resonances than higher 
mass planets if the Type I migration rate is correct. We also see 
that planets undergoing Type I migration will be migrating too fast 
to capture smaller bodies into any second-order resonances, while 
planets undergoing Type II migration will be able to capture smaller 
bodies into these resonances. 

Next, we remark that the present libration amplitude of res- 
onant bodies might have been imprinted at the instant of cap- 
ture. First we consider how to achieve high libration amplitudes. 
For low eccentricity, bo th the Hamiltonian model and full nu- 
merical simulations ( Rei n & Papaloizoi3l2009h predict small libra- 
tion amplitudes. lAdams et al. 1 20081) argue that turbulence in a 
disc, as well as breaking resonances completely, can lead to res- 
onances with large libration amplitudes of > 60°, as seen in 
H P 128311, HP 82943 and H P 73526. Similar results were found 
bv lRein & Papaloizoij ( l2009b . Here, we have shown that large li- 
bration amplitudes can also arise from high eccentricity at the in- 
stant of capture. For Jupiter-mass planets and first-order resonances 
the critical eccentricity to achieve such libration amplitudes is of 
order 0.1 (see figure |5] and equation [3j, and it is currently unclear 
whe ther a disc could pump planetary eccentricity to this level (see 
e.g.. lMoorhead & Adamsl2008l.lBitsch & Klevl2O10l) . However, for 
lower mass planets the necessary eccentricity is lower, and it may 
not be necessary to invoke turbulent fluctuations to explain a high 
libration amplitude if such low mass resonant planets are found to 
have high amplitudes of libration. 

We can also consider how to achieve very low libration am- 
plitudes. GJ 876 hosts two planets in a 2:1 resonance with both 
resonant arguments having very low libration a mplitudes of ~ 5° 
jLee & Pealell2002h . While iLee & Peald ( l2002h obtained limits on 
planetary eccentricity at the time of capture of eb < 0.06, ec < 
0.03 for planets b and c by requiring that capture into high order 
resonances not occur, we suggest that it may be possible to improve 
this: if the libration amplitude were very low at the time of capture, 
our results show that the eccentricities would have to be eb < 0.01, 
ec < 0.02. 

These results ignore any subsequent evolution of the libra- 
tion angle following capture. Libration amplitude decreases slightly 
with continued migration, but only very weakly {A.6 oc t^^^^; see 
equation |6l, so eccentricities at the moment of capture could have 
been slightly higher, depending on the extent of the migration. We 
also note that any eccentricity damping mechanism may reduce li- 



^ We restrict ourselves to classical Type I migration; the model can be 
used for the red uced or reversed migration found for example in non- 
isothermal discs I'Kley et al."2009') or isothermal discs in the non-linear 
regime jPaardekooper & Papaloizou 2009) . so long as a numerical or an- 
alytical migration rate is known. 



bration amplitudes, and we have already mentioned that stochastic 
fluctuations can increase it. 

There may be some caveats associated with using the simpli- 
fied Hamiltonian model to model capture in a hydrodynamic disc. 
However, the full hydrodynamical simulations and those using an 
N-body integratio n with a prescrib ed migration rate give compa- 
rable results (e.g., iKlev et alj|2004l) . Since we have demonstrated 
the accuracy of our model against N-body integrations, we should 
not expect many significant differences between the Hamiltonian 
model and full hydrodynamic models. We also note that a large 
planet can open a gap in a disc and hence prevent oth er bodies from 
migrating into resonances that lie in t he gap (Pier ens & Nelsonl 
l2008l : iPodlewska & Szuszkiewic3l200 9'). which will prevent cap- 
ture into those resonances. The capture process may also be af- 
fected by damping mechanisms such as gas drag or collisional 
damping, although these may be more important before and after 
resonance passage (for example, driving bodies towards the low- 
eccentricity regime prior to capture) than during capture itself. 

Finally, we point out that our results may be useful for plan- 
etary population synthesis models which require a simple alter- 
native to N-body integrations for planet-planet interactions (e.g., 
Ilda & Lin.20ia) . 



6.3 Debris disc structure 

6.3.1 Planetesimal discs 

The non-axisymmetric structure seen in several debris discs has 
been ascribed to the migration of a planet into a planetesimal disc 
(e.g., W03, Reche et al. 2008). We showed in the previous Section 
that the Hamiltonian model reproduces well the results of W03 who 
investigated this precise case. 

Resonances affect the observed disc inorphology because par- 
ticles with nonzero eccentricity trapped in a resonance spend most 
of their time at s pecific longitudes relative to the planet (W03; 
iReche et al.ll2008l) . This gives the appearance of clumps at these 
locations when the disc is imaged, with different resonances giving 
different numbers of clumps. They are seen so long as the libration 
amplitude is low enough that the clumps are not smeared out. If 
libration amplitudes are too high, the resonant dust particles appear 
to form a ring, axisymmetric save for a gap where the planet lies. 

The Hamiltonian model is particularly useful for studying the 
resonant signatures of planets in debris discs because the fate of 
tens of thousands of particles encountering any first or second or- 
der resonance can be quickly determined, regardless of migration 
rates or particle eccentricities. Hence, many images of discs result- 
ing from the migration of different sized planets into discs in differ- 
ent states of dynamical excitation can be easily made. Thus it will 
be particularly useful for exploring the disc structures generated 
by a wide range of possible perturbers for comparison to images 
expected in the near future from such projects as ALMA. 

We now return to the case of a planet at around 40 AU orbiting 
a star of 2.5 Solar masses. This is of interest because these parame- 
ters pertain to a hypothetical planet that may be imposing structure 
on the debris disc around Vega (W03). We shall discuss whether 
the planet causes a detectable signature in the disc if it migrates at 
0.5 AUMyr~^, for varying planetary mass and eccentricity (as dis- 
cussed above, we treat the planet's eccentricity as though it were the 
particles'). Figure [13] shows the effects of varying planet mass and 
eccentricity on the probability of capture and libration amplitude 
for 2:1 resonance, which have direct relevance for disc morphol- 
ogy. We see that libration amplitude increases significantly with ec- 
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Figure 13. Libration amplitudes in degrees (black lines) and capture prob- 
abilities (red lines) for planetesimals encountering the 2:1 resonance of a 
planet, orbiting a 2.5 Solar mass star, migrating outwards at 0.5 AUMyr"^ . 
Planet masses are shown in the legend. 



centricity, as expected from figure |5] The libration amplitude does 
not, however, depend strongly on planet mass. The capture proba- 
bility, on the other hand, does, with the 1 and lOM® planet cap- 
turing very few particles at any eccentricity at this migration rate, 
although capture would occur if migration were slower. Thus, we 
can predict that a 1 or lOM® planet migrating at 0.5AUMyr~^ 
will not create a detectable signature in a disc as its resonances fail 
to capture particles. Higher mass planets, however, will create sig- 
natures, with the observed structure depending on the eccentricity. 
As the eccentricity increases, capture probability goes down (Fig- 
ure |2j and libration amplitude increases (figure |5j. Thus, at low 
eccentricities, the planet will capture many particles with low libra- 
tion amplitudes, which will cause well-defined clumps to be visible 
in the disc. At higher eccentricities, the number of particles cap- 
tured will decrease and their libration amplitudes increase, leading 
to progressively weaker reso nant signatur es at higher eccentricities. 
This explains the results of iReche et al.l I2OO8), who investigated 
this scen ario with N-b ody integrations (see their figure 10). We also 
note that iReche et al.i (2008.) found that there was little difference in 
observed structure regardless of whether the planet or the disc par- 
ticles were eccentric, further justifying us modelling the planet's 
eccentricity as belonging to the particles. 

It is worth noting that converting between the observed 
dust morphology and the underlying planetesimal population is 
not a trivial problem, since the resonances can affect the colli- 
sional evolution of planetesimals ( IWvattll2006l : lOueck et al.ll2007| ; 
[stark & Kuchner 20Q9), while dust grains can be liberated from 
a resonance b y having a h igher radiation pressure coefficient than 
their parents ( IWYatlll2006i) or acquiring significant velocities rel- 
ative to th eir parents as a re sult of the collisions which gener- 
ated them jKrivov et al.ll2007t) . A fully consistent disc model will 
need to take these effects into account. Since the dynamical be- 
haviour of particles in the Hamiltonian model can be converted 
into physical positions and velocities, the Hamiltonian model can 
be used to generate a collisionless seed distribution for a collisional 
grooming algorithm to couple dynamical and collisional evolution 



jStark & Kuchnej|2009l) . Currently such models use collisionless 
N-body integrations to generate the seed distribution. In the case of 
particles experiencing Poynting-Robertson drag, the Hamiltonian 
model will enable seed distributions for many particle sizes to be 
quickly generated, allowing a finer sampling of the dust size distri- 
bution. 

These results may be applied in the Solar System to consider 
the implications of the observed libration widths of resonant KBOs 
for their eccentricities at the time of Neptune's migration. Figure|5] 
and Equations |7] and [13] can be used to assess individual objects. 
However, a general p oint we can make already is that the observed 
high libration widths i Lvkawka & Mukai 2007) would imply a high 
level of excitation, as already suggested bv lChiang et al.l ( |2003|) . 

These results can also be applied to consider the question of 
whether a disc can be stirred by a planet migrating towards the disc 
but too rapidly to capture bodies, which is known to be possible for 
diverging orbits (e.g.. , Chiang et al..2002) . 



6.3.2 Dust discs 

For less massive debris discs, such as the Solar System's zodiacal 
cloud and anticipated extrasolar analogues, particles can migrate 
substantial distances under PR drag before being destroyed by col- 
lisions. Note that currently observed debris discs are collision dom- 
inated andso the necessary migration cannot come from PR drag 
( IWvattll2005l) . Here we consider the scenario where the migration 
is due to PR drag acting on the dust, while the planet remains fixed. 

As discussed previously, dust drifting under PR drag will en- 
counter a succession of resonances. This is illustrated in FigurefTTI 
The right hand axis of ordinates shows the dust particle size corre- 
sponding to a particular migration rate, given hy D = 1.4L*/ (pa), 
with the stellar luminosity Li, measured in Solar luminosities and 
the dust density p measured in kgm^' (Wyatt etal. 19990. For 
p = 2500 kgm~^, FigurefTTI predicts that a IM® planet at 1 AU 
will be unable to capture low eccentricity particles of size less than 
about 10 micron into any resonances, but would capture particles 
of size 10-100 microns into high j resonances, and particles of size 
> 100 microns into only the 2:1 or 3:2 resonances. 

In Figure [14] we show the capture probabilities for 12 micron 
and 120 micron dust grains drifting under PR drag into the Earth's 
and Mars' resonances. Eccentricities were drawn from a uniform 
distribution between and 0.2; the particles originate in th e aster- 
oid b elt where the mean eccentricity is 0.15 ( Murra y & DermottI 
I1999I) but will have had their eccentricities excited by crossing the 
secular resonance and damped by PR drag as they migrate. Par- 
ticle sizes were drawn uniformly in the ranges [11, 13] microns 
and [110, 130] microns; note that the size distri bution of the zodi- 
acal d ust at 1 AU peaks at around 100 microns dLove & Brownie j 
I1993I) . As expected, more of the larger particles are trapped than 
of the smaller ones. For particles of around 12 microns, with this 
eccentricity distribution, the Hamiltonian model predicts that the 
most populated resonances with Earth are the 5:4 and higher j, 
with lower j resonances being less populated. This agree s qualita- 
tively with the numerical study of IPermott et al.l ( 1 19941) . who in- 
vestigated trapping of 12 micron grains originating in the asteroid 
belt into Earth's resonances with an N-body code. They found that 



* Note that because the PR migration rate depends on the particle's semi- 
major axis, a given migration rate coiTesponds to a slightly different particle 
size at different resonances. For this reason, the right-hand axis shows not 
particle size but size/o. 
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Figure 14. Capture probabilities for 12 and 120 micron dust grains migrat- 
ing into the Earth's and Mars' resonances, assuming particle eccentricities 
are randomly distiibuted in the range [0, 0.2]. 



very few (< 1%) particles are captured into the 4:3 and 5:4 reso- 
nances with Earth, rather more (2-3%) into the 6:5 to 10:9 reso- 
na nces, after which cap ture probaility decreases again (see Fig 4a 
of iDermott et al.l [l994h . Our capture probabilities are somewhat 
higher, which may be due to different eccentricity distributions or 
our neglect of stellar wind drag, but the picture of preferential trap- 
ping into j ~ 6 and higher is the same. 

Mars, being less massive than the Earth, traps fewer particles 
than the Earth. Typically we find it only captures around 10-40% 
of those Earth captures, depending on particle size and resonance j. 
This may explain why Mars lacks an observed resonant dust cloud 
while the Earth does not: too few p articles are captured to yield a 
detectable signature. iKuchner et al.l (2000) searched without suc- 
cess for a cloud trailing Mars, finding an upper limit for its frac- 
tional overdensity of 18% of Earth's cloud's overdensity, broadly 
consistent with the relative numbers of particles we find trapped 
into resonance with the two planets. If we had given particles en- 
countering Earth's resonances lower eccentricity, which may be ex- 
pected as they have had longer to be damped by PR drag. Mars 
would capture even fewer particles relative to the Earth. Martian 
capture probabilities may be further reduced due to the influence of 
the corotation resonance (W06), since the planet's eccentricity is 
the same order as that of the particles. A full investigation of this is 
beyond the scope of this paper, but it seems that Mars should have 
a resonant ring structure similar to that of Earth, at a level that may 
be close to the limits placed by current observations. 

We intend to conduct a more thorough investigation of debris 
disc structure using our Hamiltonian model in a subsequent paper. 



6.4 Limitations 

We have included only one resonant term in our integrations. 
iHolman & Murray] ( 1 19961) showed that resonances can be unsta- 
ble when the additional resonant term associated with an eccentric 
planet is important, and Q06 showed that the probability of being 
captured can be reduced by the additional resonant term. However, 



since iReche et al. I j2008l) found similar results in full N-body in- 
tegrations when using eccentric particles and eccentric planets, in- 
cluding the extra resonant term explicitly may not be necessary. 

We have restricted our numerical integrations to the case of a 
single planet plus a test particle. This may not adequately model 
two-planet systems where the planets are of comparable masses. In 
this case the equations given in the Appendix. |A16| and |A28l can be 
used. However, the increased number of free parameters will make 
a numerical study with a comparable grid resolution less practical. 

We have neglected any eccentricity damping, which can oc- 
cur together with change of semi-major axis when migrating in a 
gas disc or under PR drag, or experiencing dynamical friction in 
a planetesimal disc. This will limit the eccentricity growth due to 
continued migration in resonance. The effect on lib ration behavio ur 
is dependent on the precise nature of the damping l lGomesiri995[) . 

We have only taken the lowest order term in the disturbing 
function, which is justified for low eccentricities. However, when 
eccentricities are high during migration, or are driven higher as mi- 
gration continues after capture, higher order terms may become 
import ant. This is the cas e for the resonance between GJ876b 
and c ( lLee&Pealell2002l) and for the two libration centres of 
the 2:1 resonance, seen for eccentricities above ~ 0.03 (W03; 
iMurrav- Clav & Chian^l2005l) . While post-capture eccentricity in- 
crease will not significantly affect the capture probability (a particle 
must cross a similar separatrix to enter the libration region in the 
first place) it will change the behaviour of libration amplitude and 
centre following capture. 

We have not considered higher order resonances such as 
the 4:1, which can be important for very massive planets 
jKuchner & Holmaij|2003h . 



7 SUMMARY AND CONCLUSIONS 

We have systematically investigated the Hamiltonian model for 
capture of a test particle into first and second order resonances with 
a planet. The model reduces the full complexity of the restricted 
three body problem to only one degree of freedom, the resonant 
angle 6, and two parameters, proportional to the particle's eccen- 
tricity and the migration rate. Only the relative migration rate of 
the two bodies, in the form di/ai — 0,2 /a2, is important. External 
and internal resonances behave the same way, but with different 
proportionality constants. 

We confirmed previous work showing that capture into reso- 
nance is certain at low (rescaled) eccentricities and migration rates, 
possible with decreasing probability at high eccentricities, and im- 
possible at low eccentricities and high migration rates. As eccen- 
tricity increases, the transition from certain capture at low migra- 
tion rate to impossible capture at high migration rate broadens. At 
higher eccentricities, capture is possible with faster migration than 
at lower eccentricities. This effect is more pronounced for second 
order resonances. 

We have also found the libration amplitudes and centres of the 
resulting resonant motion. The libration centres are offset from 6 — 
TT, the centre in the absence of migration, by an amount increasing 
with the migration rate, agreeing with previous work. In addition, 
we have found that the offset is almost independent of eccentricity 
for first order resonances, except at high eccentricities. Libration 
amplitudes are small if capture occurs at low eccentricity and low 
migration rate. They are somewhat larger for migration rates just 
less than critical. Very large libration amplitudes (> 90°) can be 
attained if the initial eccentricity was very high. 
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We have also found the jumps in eccentricity when a particle 
encounters a resonance but is not captured. In the case of slow mi- 
gration, the eccentricity is always driven down, in agreement with 
adiabatic theory. When the migration is fast and the eccentricity 
low, eccentricity jumps up. However, when the migration is fast 
and eccentricity high, eccentricity can jump either up or down. 

We c hecked the res ults of the model against the N-body simu- 
lations of Wyatt' ('2003'), finding excellent quantitative agreement 
for capture probabilities, and qualitative agreement for libration 
amplitudes and offsets. We found that accounting for the particles' 
eccentricity is necessary to fully explain the dependence of capture 
probability on migration rate. While the effect of eccentricity has 
long been understood in the context of the Solar System, its im- 
plications for extrasolar planets are less well explored, and for low 
mass planets in particular its effects are important. 

We have applied our model to several situations in which 
planet or particle migration is likely to occur. We can easily de- 
termine whether capture occurs for planets migrating in Type I or 
Type II regimes. The pre-capture eccentricity can be constrained 
by the present libration amplitude, with higher eccentricities giving 
higher libration amplitudes. For planets migrating through a gas 
disc, a non-zero eccentricity prior to capture can lead to large li- 
bration amplitudes such as those seen in the HD 128311 system. 
We find that, if planetary eccentricity can be raised to e > 0.1 
for Jupiter-mass planets, for example by planet-disc interaction, the 
resonant capture process by itself will result in high libration am- 
plitudes without the need to invoke extra mechanisms such as tur- 
bulent torques. For lower mass planets the necessary eccentricity 
is lower. Also, very low libration amplitudes (e.g., GJ 876) sug- 
gest a low eccentricity during the capture process. We also predict, 
based on classical formulae for migration rates, that lower mass 
planets will be found in higher j first order resonances than higher 
mass planets, and planets migrating in the type I regime will be 
moving too fast to capture smaller particles into any second-order 
resonance. The model may also be useful for population synthesis 
of multi-planet systems where it is necessary to account for planet- 
planet interactions without recourse to full N-body integrations. 

We then discussed debris disc structure. A planet that has mi- 
grated into a disc will impose clumpy structure on a dynamically 
cold disc, but the clumps will be at a lower level and smeared out 
for migration into excited discs. For dust migrating under PR drag, 
the model can explain the structure of the Earth's resonant ring rea- 
sonably well. We predict that Mars has a dust ring at a level < 25% 
that of Earth's, consistent with observed upper limits. 

The data on capture probabilities, libration amplitudes 
and offsets, and eccentricity jumps, are available on-line at 
[htt p : / / www .ast.cam.ac.uk/~ajm2 33/ and at the jour- 
nal website and may be used provided that this work is cited. 
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of the Hamiltonian model such as the evolution of libration ampli- 
tude and offset with time. Although many of these results may be 
found in the literature, we think it useful to gather them together 
for reference. 

First we derive the scalings to link physical variables such as 
eccentricity to the dimensionless ones such as the momentum J. 
This closely follows Q06; we generalise to consider two massive 
planets. 

We use Poincare's canonical variables; the generalised mo- 
menta are 



(Al) 



which are conjugate to the mean longitudes A; and the negative of 
the longitudes of pericentre vji respectively. Here, yLi = Q {nit, + 
rrii). The Keplerian part of the Hamiltonian is therefore 
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We work with low eccentricities, sothatFi/Ai ~ e^/2. 

We now non-dimensionalise as follows. Distances are put in 
units of ao, the semi-major axis of the inner planet when the nom- 
inal resonance is reached. Times are put in units of the inverse of 
the mean motion of the inner planet at that point. Masses are put in 
units of the mass of the inner planet. With this choice of units, the 
Keplerian Hamiltonian is 
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with primes denoting the non-dimensional quantities. We now drop 
the primes, and make the approximation /i2 ~ Mi = 1, so that 

1 ..3 
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where v — m2 /mi is the ratio of the two planets' masses. 

To study the dynamics as the system passes through reso- 
nance, we expand about the nominal resonance location where 
Ai = Ai,o, with Ai,o = 1 and A2,o = ua^^^'^, qq = [{j — 
being the semi-major axis ratio at the nominal resonance. 
This gives the Keplerian Hamiltonian 

ttkep = /i - -h H : — h - - — : — I u h 

2 J 2 V J 
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(A5) 



where li = Ai ~ Ai,o measures the distance from resonance. This 
expansion assumes l2a^^^ /v <C 1. This is satisfied for all f since 
I2 — 0[v) as — >■ 0. We just require the bodies to be close to 
resonance in semi-major axis. 

From the disturbing function we retain the relevant resonant 
terms. The resonant Hamiltonian is 
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(A6) 



APPENDIX A: HAMILTONIAN MODEL 

Here we detail the mathematical derivation of the scalings from 
physical variables to dimensionless variables, and some properties 



Again, the r coefficients are functions of a and regarded as con- 
stant. In terms of the fi functions, they are 
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Table Al. Values of a and fi for first-order resonances. 

where 1(2 : 1) = 1 for the 2:1 resonance and is zero otherwise; 
this is from the indirect disturbing function. 

The fi functions are tabulated for first-order resonances in Ta- 
ble i IaB . 

Al First order resonances 

Now we consider first-order resonances. There are two resonant 
angles 9i and 6*2; we effect a point transformation to these angles: 

Ql = jA2 + (1 — — 7J7l 

d-Z = jA2 + (1 — i)Al — W2 
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The new momenta satisfy 

(l-j)(Jl+ J2) + J3 = 
j(Jl + J2) + J4, = 
Jl = 
J2 = 

The Keplerian Hamiltonian is then 

■Hkcp = -3(1-jV3(Ji+'^2) 

-^J4(Jl + J2)- 
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The resonant Hamiltonian is 

"Hrcs = rifiJ2^^ 008 6*2 + J"!,! COS 6*1 . 

This reduces the system to two degrees of freedom: the mean lon- 
gitudes 9^ and 64 do not appear in the equations. The momenta J3 
and J4 are therefore constants; we dropped terms involving only 
these momenta from Equation JAlOt . We simulate migration by 
explicitly varying the momenta J3 and J4 with time. We have 
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Al.l Reducing number of parameters 

We return to the general Hamiltonian (Equations I A lOtjA lit . This 
may be written 

H = aJi + PJ2 + 7 J? + 27J1 J2 + 7JI 
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These six parameters we reduce to four by rescaling momentum: 
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Note that 7 < 0, and ri.i > for first-order resonances. 

We may write the parameters in terms of physical quantities 
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3(j - 1) (\/a2/a2,o - + "^mV^ h 



2,211/3 



31/3 (w/27f/' [(1-JT + ^] 

_a^^^i3^^2a^^^I(2^ 



,2 -2 1 1/3 



t 



(A19) 
(A20) 



Migration entails varying the coefficients a' and /3' with time. 
We have: 



(A12) ^ ^ ^ 



5I/3 / 02 
\ 02 



/ / \ 4/3 / \ 2/3 

20-- 1)1/3 (1 + 4^) 



.(A21) 



A7.2 Limiting case 1: outer test particle and zero eccentricity 
planet 

We consider the limit v ^ Q. This is relevant for either a planet 
capturing planetesimals as it migrates outwards, or dust spiralling 
inwards under P. R. drag. In this case a more convenient scaling is 
using ri.o rather than ri,i. The non-constant terms in the Hamilto- 
nian are 
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where 



3(j - 1) (^^02/02,0 - ^/af^ + 2^.al'^ f2 



(A22) 



(A23) 



Note that all terms in Equation | |A22I ( are 0(1/"). For a migrating 
planet, the migration rate is given by 

d/3o ^ 3'/^(j-l)(a2/a2,o-ai) 



(A24) 









fro 
/53 


3 


0.480750 


0.598759 


1.98591 


5 


0.711379 


3.27381 


5.68728 


7 


0.799064 


7.87052 


11.3173 


9 


0.845740 


14.3866 


18.8674 


11 


0.874782 


22.8216 


28.3365 


13 


0.894608 


33.1755 


39.7246 


15 


0.909009 


45.4473 


53.0314 


17 


0.919944 


59.6382 


68.2557 


19 


0.928532 


75.7482 


85.3994 


21 


0.935455 


93.7765 


104.462 



Table A2. Values of a and fi for second-order resonances. 



A1.3 Limiting case 2: inner test particle and zero eccentricity 
planet 

The limit v — ^ 00, — > 0, pz^ = m2/mo constant, is relevant 
for a planet capturing planetesimals as it migrates inwards. The 
Hamiltonian is 



(A25) 



(A26) 



n — Ji + Qoo Jl — Jl ' COS Oi 

where 

3(j — 1) (^V 02/02,0 — \/aT) + 2/i!^ao/2 

3i/3[/,,.ao/27(l-j)]'''' ■ 

All terms in Equation l lA25t are O The migration rate in this 
case is given by 

daoo _ 3^^^(02/02,0 - ffli) 



(All) 



A2 Second order resonances 

For second order resonances there are three resonant terms. The 
scaled Hamiltonian (cf. Equation lA16t is 

n' = a'Ji + p'j^ + {j[ + j:,f 

+s'Ji^^J2^^ cos(6li - 62) - Jl cos 201 
+r[Ji^^J2^^ cos(6»i - 62) + ra J2 cos26l2, (A28) 
with the rescaled time 

t' = 2fiuaof45t, (A29) 
the rescaled momenta 

3[(2-jf + 4^] 



and the parameters 



^2 = 



3 [(2 - j)ay^ + jag^/'] + 4^^/00/2 
A fiuaof 45 

3 [(2 - i)a'/' + jaUl^^] + A^al^\f2 
4 fii/aof 45 

1/4 r 

«o 749 

[2/53 - f aoI(3 : 1)] 



21^/45 

The migration rate is given by 



(A30) 

(A31) 

(A32) 
(A33) 
(A34) 



da' 
dF 



d^' 

dF 



3(j-2)(d2/a2 



- ai 



(A35) 



The values of the ji for second-order resonances are tabulated 
in Table | |A2] |. 

For an outer test particle and a planet on a circular orbit we 
instead rescale by r2,o, giving 



t' = 



3/2 
«o M 



27 



2,f53 - — aoI(3 : 1) 



d^ 
dt' 



3i' 



16m [2/53 - f aoI(3 : 1) 
3 



2 

-62 



2q;^''V [2/53 - f aoI(3 : 1)] 
3(j - 2)(a2/a2 - ai) 
4alp? [2/53-faoI(3: I)]'' 



(A36) 
(A37) 

(A38) 
(A39) 



A3 Evolution in resonance 

A3. 1 First-order resonances 

After a particle has been captured into resonance, the trajectory 
continues to evolve so long as the imposed migration is still present. 
An analytical description of the behaviour is possible when |/3| is 
sufficiently large. Both the amplitude and the offset of the libration 
centre of the libration in 6 decrease with time, while the momentum 
J increases. 

For first order resonances, the libration centre can be found by 
considering the equations of motion for 6 and J: 



e = 2J + /3- -J^^/^cose 

j = J^'^sinO, 

whence, 

e = 2J^/^sin6' + /?+ ij"^sin6'cosf 



(A40) 
(A41) 



1/2 



2J + /3--J- 



-1/2, 



cos 6* sin 9. 



(A42) 



Now, the elliptical fixed point occurs at approximately Jcq = 
—13/2 for large thus giving the increase in momentum. For 
large |/3| and small 9, Equation lA42l tells us that 



so that the libration centre offset is given by 



(A43) 



Resonance capture 1 9 



(A44) 



To find the evolution of the libration amphtude, we consider 
the equations of motion in Cartesian coordinates, following the 
canonical transformation 



X — \/2J cos ( 



y 



2 J sine. 



(A45) 



We calculate the libration amplitude by imposing a constant area 
on the small elliptical trajectory enclosed by the orbit; a knowledge 
of the aspect ratio of the ellipse then suffices to determine the am- 
plitude in y, and hence 6. 

We improve the estimate for the equilibrium point by a single 
iteration of the Newton-Raphson methodic giving 



Xeq~a + l/(2y2a2), 
where a = \/—l3. Writing X 
X = -y/V2a, y = 2a^ 



(A46) 

a;cq and linearising gives 

(A47) 



to leading order in a. This describes motion in an ellipse about 
= 0, with ratio of semi-axes AX/ Ay = 2~^^'''a~'^^^ . The area 
enclosed by a trajectory is thus A oc a~^^^{Ay)^. Since A6 ~ 
Ay/x, and a; ~ a, we find 



A6I ~ Ca 



-1/4 



OC P' 



-1/8 



(A48) 



for some constant C. This constant is not determined by this anal- 
ysis; we obtain it numerically in §|4] 

This also gives the period of libration: 



-1/4 



(A49) 



A3. 2 Second-order resonances 

For second-order resonances, a similar analysis gives 

- ^ 

^''^ 4/3' 

Ae^m-'^\ 

and 



(A50) 
(A51) 

(A52) 
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° Note that working directly with Joq : 
the equation for X vanishing. 



-f}/2 leads to the linear terms in 



